Origin of second- harmonic generation in the incommensurate phase of K2Se04 
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We show that a ferroelectric phase transition takes place in the incommensurate phase of the 
K2Se04 crystal. The ferroelectric character of the IC phase explains the second-harmonic generation 
observed in the corresponding temperature range. 

PACS numbers: 64.70.Rh, 64.60.-i 



The crystals of the A2BX4 family (about 20 crystals )i 
undergo two successive phase transitions on cooling: an 
incommensurate (IC) transition at temperature Ti with 
a modulation vector near the &/3 point in the Brillouin 
zone, and upon further cooling at temperature Tc a com- 
mensurate lock-in transition to the triple-period ferro- 
electric phase with the modulation vector 6/3. The most 
striking peculiarity of the K2Se04 crystal, which belongs 
to the A2BX4 family, is the second-harmonic generation 
(SHG) in the IC phased of the same intensity magnitude 
as that in the 6/3- phase. Low temperature 6/3- phases 
of the A2BX4 type crystals are improper ferroelectrics, 
and therefore the SHG is a normal manifestation of this 
feature. However, so far there was neither pointed out a 
plausible reason for the SHG in the IC phase of K2Se04, 
nor explained why it takes place not in all the crystals of 
the A2BX4 family. It is worth noting that the SHG was 
observed also in the IC phase of quartz by Dolino and 
Bachheimer (1977)- and in the IC phase of the ammo- 
nium fluoroberyllate (NH4)2BeF4 by Alexandrov et al. 
(1978)^, and both observations also can be explained us- 
ing the theory developed below. The attempt to explain 
the SHG in the IC phase of (NH4)2BeF4 was carried out 
by Golovko and Levanyuk^ based on the spatial disper- 
sion of the dielectric constant. However, the expected 
effect appeared to be small, and besides, it is not ex- 
plained why in such a case the SHG is not observed in 
all the crystals with the IC phases. 

In the present paper we show that within the IC phase 
of the K2Se04 crystal, necessarily a transition to the fer- 
roelectric IC phase takes place. We analyze the exact 
solution in frame of the Landau theory. The physical 
reason for such a ferroelectric transition is the coupling 
between the crystal's polar symmetry vibrational mode 
P and some displacements of the IC domain walls. We 
demonstrate that relative displacements of the domain 
walls in the domainlike IC structure induce polarization 
along the z direction, and can be viewed as an additional 
phonon branch of a polar symmetry. This mode is not 
explicitly introduced in the equations below. It is a com- 
ponent of the IC phase 'p{x) function, and is taken into 
account implicitly, in the integral form. The coupling of 
the two modes renormalizes the frequency of the lower 
mode down to zero within the IC range, which induces a 
ferroelectric transition. 

Structure of the IC modulation of K2Se04 was defined 
in the neutron diffraction by lizumi et al& The thermo- 



dynamic potential describing the IC phase of K2Se04 
was detailed studied by Sannikov and Levanyuki For 
the two-component order parameter with components 
Vi — ^70 cosM and 772 = ijq sinw it can be written as: 
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Bvo + Fvl + "^Pvl cos 3w + ^-p-^ + hP^ - P£} (1) 

where P is polarization vector in the z direction, £ is 
electric field in the z direction. In case of cr = 0, Z? > 
this potential would describe a phase transition from the 
symmetric high-temperature phase to the triple-period 
commensurate phase at the temperature a — Q. For 
/ < the equilibrium values of the modulation phase 
u are 0, 7r/3, 27r/3, . . . (otherwise 7r/6, 37r/6, . . .). How- 
ever, in virtue of the Lifshitz invariant (ct 7^ 0) the phase 
transition takes place from the high-temperature phase 
to the IC phase at a > with u = kx + ip and k = —a/D. 

Substituting u ~ kx+(p{R) in potential (1) one obtains 
for the if- and P- dependent part of potential (1) the 
following expression: 

$ = J dx{jvU^f + {vocosi6kx + 6ip) + 

^Pr]^ cos{3kx + 3ip) + ^-P^ + hP* - P£} (2) 

Minimization of this expression with respect to ip gives 
the equation: 

D-—r + 3/ryn sin(6fca:: + (Sip) H F770 sin(3fca; -|- ^ip) = 

ox'^ 2 

After single integration over x and taking for the inte- 
gration constant the form |/|?7o(l + c) one obtains: 



^ = Vl/l7D'7oWl + c-cos6u+|y|-3Cos3u, (3) 

where u — kx + ip. The integration constant c, as well as 
the polarization vector P, can be defined via minimiza- 
tion of the potential with respect to these parameters 
after substituting the solution of Eq. (3) in the poten- 
tial (2). Prior to considering the ferroelectric transition 
in this system let us discuss the solutions of Eq. (3) for 



2 



12 
10 
8 
6 
4 
2 



kx+(pix) 



a) 



II 



III 



10 20 30 40 50 60 70 
kx+(p(x) ^ / 




10 20 30 40 50 60 70 

FIG. 1: a) Coordinate dependence of the IC modulation phase 
for different temperatures. The curves I , II, III represent 
various stages of evolving of the domainlike modulation on 
cooling: curve I, j^q ~ 0.21, c — 0.138; II, rjQ = 0.23, c — 
0.010095; and III, ryg = 0.2326288, c = 9 x 10"^ (the last 
is very close in i]o to the lock- in phase amplitude rjo). The 
coordinate x is given in multiples of the lattice parameter. 
The height of each step on the curves is equal to 7r/3. Each 
domain (plateau) has nonzero polarization vector in the z di- 
rection. However, the neighboring plateaus are polarized in 
opposite directions, providing crystal's zero integral polariza- 
tion, b) Transformation of the curves I , II, III induced by 
the electric field, applied in the z direction. The effect of the 
electric field £ comes to the increasing of the plateaus' areas 
if the polarization vector is in the direction of £, and to the 
reduction of those otherwise. The height of each step remains 
7r/3. Such a transformation should take place spontaneously, 
without an external field. 



the case of P = 0. The solution to this kind of equa- 
tions for the IC structures was analyzed in the pioneer 
work by Dzyaloshinski^ (it was also detailed studied by 
Golovkc^), and for Eq. (3), it has the form: 



ip{x) + kx 



iam(3/[fi7D,S, 



c 



(4) 



where am is the elliptic amplitude function and k is 
the IC vector at the temperature of the IC transition. 



(In the present paper the notations of the Wolfram Re- 
search, " Mathematica" program for the eUiptic functions 
are used.) Substituting solution (4) in the potential (2) 
and minimizing over the parameter c one obtains: 



Trfc 



(5a) 



E is the complete elliptic integral of the second kind, 
which has the asymptotic 



V^E(— ) 
c 



2(1 --Inc) 



for c — > 0. 



Eq. (5a) can be numerically solved with respect to c for 
given coefficients /, D, fc, t/q, and one obtains that the 
parameter c monotonously decreases from +oo at rjo = 
(i.e., at the IC transition Ti) down to c = at the 
continuous lock- in transition Tc- As it follows from Eq. 
(5a), at the lock-in transition point the IC amplitude is 
tending to 



V27rfc 



VD/\f\ with c ^ 0, 



(5b) 



which gives that k ^ rj^ , if the remaining coefficients are 
of the same order of magnitude. 

In Fig. 1(a) the solutions given by Eq. (4) for k — 
27r/30, D/\f \ = 1 and P for the three values of the IC 
amplitude are plotted. At the IC transition Ti (c — > oo) 
the corresponding graph must look like a straight line 
with slope k [if = const, i.e., sinusoidal modulation), and 
on cooling it acquires the domain-like features. Due to 
the symmetry, each point x on the curves of Fig. 1 (a) is 
associated with a nonzero polarization, which is spatially 
distributed as P ~ 77g cos(3A:a; -I- 3(^) with zero-integral 
effect. Meanwhile, each domain there has a nonzero- 
integral polarization. In the areas of any of two neigh- 
boring domains -\-kx n7r/3 and ip -\- kx ~ (n + l)7r/3 
respectively, giving rise to the P-s of opposite signs. 

In Fig. 1(b) the transformation of the curves of Fig. 
1(a) under the applied electric field £ is shown. The field 
£ is taken into account via replacing in Eq. (3) of the 
polarization vector P by P = Xo£- An identical trans- 
formation of the domains takes place in the presence of 
nonzero polarization P in crystal (when £ — Q). Corre- 
sponding solutions were found both analytically and nu- 
merically (using the "Mathematica" software) with the 
same result. Note that the integral polarization vector 
contributed by the nonequal domains in Fig. 1(b) is not 
zero. 

Now we show, that the above described transformation 
of the domain structure from that shown in Fig. 1(a) to 
Fig. 1(b) takes place spontaneously, i.e., without apply- 
ing an external electric field. In other words, there takes 
place a ferroelectric phase transition inside the IC phase. 
Let us calculate the dielectric susceptibility of the crystal 
in the IC phase: 
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Differentiation of potential (2) gives: 
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3P2 



p=o 



_ 1 3rr/^ 



dx, (6) 



p=o 



t = 27r/g is the IC period, given by Eq. (10) below. In or- 
der to calculate the derivative {du/dP)p^Q, we introduce 
notation 



Hu, p) 



^1 + c — COS &v + ijj^ cos"3w 



Integration of Eq. (3) gives in this notation, F(u, P) — 
^J\f\/D'qlx = 0, from which follows that d¥/dP + 
{d¥ /du) du/dP = 0. The latter allows one to calculate 
the derivative {du/dP)p^^ via calculation of the corre- 
sponding derivatives of the function F in the integral form 
(see the Appendix). After substituting of the calculated 
derivative in Eq. (6), one obtains: 



^2$ 

ap2 



P=0 



Xo - 



Xo - 



27r/3 



sin ^udu 



At\f\ir}QC J \/l + c- cosSm 



Dr^ 



4ij?2|/|2C2 



2Fi|i,|,2, (7) 



where 2F1 (1/2 , 3/2 , 2, — 2/c) is a hypergeometric func- 
tion. This hypergeometric function has the asymptotic 



2 ' 2 



for c — > 0, 



which well fits the function in all the range < c < 1, 
i.e., not only at c ^ 0. As one can see from Eq. (10), the 
IC period has an asymptotic: c 0, t ~ |lnc|/?7o 
shown below, on approaching the lock-in temperature Tc 
the coefficient c is tending to zero as c| lnc| ~ \T — Tc\). 
So, for the susceptibility x one obtains: 



X 



Xo 
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16|/|cK 



Xo 



|/|c|lnc| 



(8) 



and the stability coefficient necessarily vanishes 

(since the c -dependent term in Eq. (8) is negative and 
diverging with c ^ 0), inducing ferroelectric phase tran- 
sition at the temperature = 0. A loss of stability 
necessarily takes place even for the case of large (i-<3-j 
even for a not soft polar mode of the high temperature 
phase), and it is expected in the medium temperature 
range of the IC phase (c < 0.1). 

However, some alternative situations are possible near 
the loss of stability x~^ = 0. As it follows from Eq. (3), 
the polarization P which appears below the transition 
should be sufficiently small: 



P<cr,l\f\/r 



(9) 



Otherwise, the expression under the root in Eq. (3) takes 
negative values for some u, what means that the crystal 
does not fall within the IC phase, but within the commen- 
surate lock-in phase instead, skipping the ferroelectric IC 
phase. In terms of Fig. 1(b), it is equivalent to the in- 
finite increasing of the one domain's size, which covers 
all of the crystal area. The right hand side of Eq. (9) is 
very small, since rjQ in Eq. (9) is a small parameter, and 
c tends to zero in the low-temperature range of the IC 
phase (see also below). In the case of a first-order ferro- 
electric phase transition giving rise to a large P (in Lan- 
dau theory, it means a negative fourth-order term in the 
potential's P- expansion), inequality (9) does not hold, 
and the transition transfers the crystal directly to the 
lock-in phase, skipping the ferroelectric IC phase. Even 
for some order-disorder transitions, which are very close 
to the second-order type (the polarization vector P below 
the order-disorder transitions is larger, due to the small 
Curie constant compared to that for the displacement 
type) the expected ferroelectric transition is more likely 
to turn into the transition to the lock-in phase, skipping 
the ferroelectric IC phase. Note that among the crystals 
of the A2BX4 family only K2Se04 is of a displacement 
type, and therefore for this crystal the ferroelectric phase 
transition (from the IC to the ferroelectric IC phase) can 
be considered as a most likely. In other words, the IC 
structure of K2Se04 becomes ferroelectric with the mod- 
ulation, spatially distributed hke that in Fig. 1(b). One 
should expect that the ferroelectric transition takes place 
in the beginning of the domainlike modulation formation, 
i.e., at temperatures near the curve I in Fig. 1(a), when 
the c coefficient falls below unity. 

As it follows from Eqs.(5a, 5b), for estimates of in- 
equality (9), the asymptotic of the parameter c — > can 
be used: 



cine « 8 



near the lock- in point we presented % — ??o ~ ^C^c — T), 
where a is some constant. So, Eq. (9) can be presented 

as P < 16a{T - T,,) T]'^\.f\/r, when T ^ (we ne- 
glected a logarithmically diverging term In c, which does 
not change order of magnitude of the estimated values, 
since, as discussed below, \Tc — T\ can not be sufficiently 
small for it). 

For obtaining of the inharmonic fourth order term in 
the P expansion of the potential one should calculate the 

derivative ^9"'$/9P^^ . We shall introduce here only 

the result of this calculation. In addition to the hP^ 
term in Eq. (2), one obtains a negative, strongly diverg- 
ing as with c ^ term. If the sum of the inhar- 
monic P^ terms becomes negative earlier than x~^ = 0, 

then takes place a first-order transition. The latter makes 
more probable the first order transition, even if the coef- 
ficient /i > is not small. 

Now we shall discuss some features of the IC phase 
for P = case in order to compare it with behavior 
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of the ferroelectric IC phase. To study the temperature 
evolving of the ferroelectric IC structure upon cooling 
one should minimize potential (2) also over the IC am- 
plitude r]o- However, we imply in the present paper only 
increasing of the IC amplitude t]o upon cooling, and do 
not discuss the character of the lock-in transition. 

First we discuss some structure peculiarities of the non- 
ferroelectric IC phase of K2Se04. The length of each 
plateau in Fig. 1(a) is six times shorter than the IC wave- 
length 2'K/q, where q is the IC wavevector [q is equal to 
the average slope of the given curve in Fig. 1(a)]. Prom 
this geometrical fact it follows, that the IC modulation 
argument u = kx + ip can be presented as u = qx + 0{x), 
where 9{x) is a periodic function with period 2'K/Qq. As it 
follows from Eq. (4), roughly 6{x) can be approximated 
by sm&qx. The IC satellite reflections observed in the 
X-ray diffraction are induced by the Fourier-components 
of the modulation function % cos[(| + q)x + 0{x)\. The 
IC satellites observed for such a modulation should ap- 
pear in the points | + f + 9 ^ 6(7, ±12(j, . . . near the 
main Bragg- reflections G. In other words, the mani- 
festation of the domain-like IC structure in K2Se04 in 
diffraction should be only the simultaneous observation 
of the IC satellites of about the same intensity (asymmet- 
rically disposed with respect to the | point) in the points 
I — 5(7 and | + 7(7 near G. Intensive IC satellites in the 
points I -t- g ± 2^, ±3(7, ±4^, iSg, which are observed 
in diffraction, are not contributed by the domain features 
of the modulation, and can not be viewed as a manifes- 
tation of a domainlike IC structure. The satellites | — 5g 
and I + 7g were not observed for K2Se04, which proves 
that the domainlike structure is not well developed in 
K2Se04, and the ferroelectric phase exists only in the 
medium range of the IC phase. So, as it follows from the 
exact solution of Eq. (3), any crystal of the A2BX4 type 
with well developed IC domain structure (i.e., sufficiently 
small c) should be ferroelectric. The ferroelectric tran- 
sition does not change the translation symmetry of the 
IC phase, and therefore the satellites disposition for the 
ferroelectric IC phase should not essentially differ from 
that for the P = case. 

We also discuss the character of increasing of the IC 
modulation period on cooling, which is logarithmically 
diverging, as it follows from Eq. (4). The IC period for 
the modulation given by Eq. (4) is 
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\f\vWc 



K(--), 



(10) 



where K(— 2/c) is the complete elliptic integral of the first 
kind. With c ^ the function K{—2/c)/y/c has the 
asymptotic 



and 



4K(--)/v^~ V2(51n2-lnc), 



VD\\nc\ |ln(r-Tc)| | ln(r - T^; 



which is logarithmically diverging. The temperature T 
can approach the lock-in temperature not closer than 
the temperature fluctuation -\/<~(AT)2~> , i.e., 

rp2 



\T - Tel > \/< (Ar)2 > and < (Ar)^ >= 



where V is the crystal volume and Cy is the heat capacity. 
So, for any crystal of a finite size {V ~ lO^^a'^), the period 
t can increase only about 10 times compared to 211 /k. in 
virtue of the logarithmical character of the divergence (it 
is also confirmed by our numerical simulations). Besides, 
in the real experiment the temperature T can not ap- 
proach the lock-in temperature Tc even closer than the 
instrumental temperature resolution AT ~ 10~^'~^/r, 
which does not allow to observe any diverging IC period, 
and the lock-in transition always should be observed as a 
discontinuous phase transition with jump in the wavevec- 
tor q down to g = 0. 

For calculation of the ferroelectric IC modulation's pe- 
riod [the .T-length of the six neighboring domains in Fig. 
1(b)], one should integrate Eq. (3) from u = through 
u = 277, which gives: 




1/1%' 



where F is the elliptic integral, and Ac = c — rP/{\f\r]Q), 
i.e., the condition Ac > coincides with Eq. (9), and 
therefore it always holds in the IC phase. The param- 
eter Ac decreases down to zero with increasing of the 
polarization vector P on cooling. Integral (11) was cal- 
culated exactly, and after being simplified to the form of 
Eq. (12) for Ac < c < 1. The first term in Eq. (12), 
which gives the size of the larger domains in Fig. 1(b), is 
diverging with Ac ^ as ~ | In Ac|. Meanwhile, the size 
of the smaller domains is given by the second term in Eq. 
(12), and it is always finite, though it also increases with 
cooling. 

Summarizing, we showed that in the IC phase of the 
K2Se04 crystal takes place ferroelectric phase transition, 
which explains the SHG in the IC phase of this crystal. 
Formation of the ferroelectric IC phase is most likely near 
the displacement-type phase transitions, though it can 
occur also in those of the order-disorder type. The IC 
period increases on cooling towards the lock-in tempera- 
ture Tc as I ln(T — Tc)|, though the size of the domains, 
polarized opposite to the prevailing polarization, remains 
finite. 



APPENDIX: CALCULATION OF {du/dP)p^Q 

In order to calculate the derivative {du/ dP) one 
should calculate the corresponding derivatives of the 
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function ¥{u, P). So one obtains 

dF 
du 



1 



p=o 



Vl + c — cos6u 



and 



that {dc/dP)p^Q = 0. Similarly {dr]Q/dP)p^Q is also 
equal to zero. So, one obtains: 



du \ _ I Off \ I I aif 



dP 
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dP 
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dl 



sin 3u 



d¥ 
dP 



2i/hgy (1 



cos 3u dv 



i6w)3/2 



sinSu 



6c|/|77|^ VI + c — cos 6m 



For integration over x in Eq. (6), the following relations 
are also used: 



ox 



27r/3 



du 



du 



3Vd 

\f\7]l J VI + c - cos6u' 



We took into account that for the equilibrium value of where we used the right hand side of Eq. (3) for idu/dx) 
c, the condition c(P) = c{—P) holds, from which follows at P = 0. 
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